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We construct exact functional renormalization group (RG) flow equations for non-relativistic 
fermions in arbitrary dimensions, taking into account not only mode elimination but also the rescal- 
ing of the momenta, frequencies and the fermionic fields. The complete RG flow of all relevant, 
marginal and irrelevant couplings can be described by a system of coupled flow equations for the 
irreducible n-point vertices. Introducing suitable dimensionless variables, we obtain flow equations 
for generalized scaling functions which are continuous functions of the flow parameter, even if we 
consider quantities which are dominated by momenta close to the Fermi surface, such as the density- 
density correlation function at long wavelengths. We also show how the problem of constructing 
the renormalized Fermi surface can be reduced to the problem of finding the RG fixed point of 
the irreducible two-point vertex at vanishing momentum and frequency. We argue that only if the 
degrees of freedom are properly rescaled it is possible to reach scale-invariant non-Fermi liquid fixed 
points within a truncation of the exact RG fiow equations. 

PACS numbers: 71.10-w, 71.10.Hf 



I. INTRODUCTION 

Recently several authors |^-]|] have used exact functional renormalization group (RG) methods to gain a deeper 
understanding of strongly correlated non-relativistic fermions in reduced dimensions. The exact functional RG yields 
an infinite hierarchy of coupled differential equations describing the change of the correlation functions due to the 
elimination and the rescaling of the degrees of freedom. For classical statistical mechanics problems the exact functional 
RG has been developed long time ago in a pioneering work by Wegner and Houghton H, who performed the usual 
three RG-steps [fflH] to derive their exact flow equation: 



1. Integrate out fields with momenta k in a shell Ao(l — dt) < |k| < Ag, where dt is infinitesimal and Aq is some 
ultraviolet cutoff. 

2. Rescale the remaining momenta by a factor 1 + dt and express all quantities in terms of the rescaled momenta. 

3. Rescale the remaining fields by a factor of 1 -f ^(1 — 77)0?^, where the anomalous dimension 7/ has to be chosen 
such that the RG has a fixed point 

Over the years several alternative formulations of the exact functional RG in field theory and statistical physics 
have been proposed p^|-p^. Note, however, that recent applications of the exact RG for two-dimensional fermions 
take only the mode elimination step 1 into account, and do not perform the above rescaling steps 2 and 3. While 
such a procedure is legitimate if one is interested in solving a given many-body problem iteratively, the Wilsonian 
RG amounts to more than that : by combining the mode elimination with a suitable rescaling of the degrees of 
freedom, the Wilsonian RG generates a mapping of the original many-body problem onto a continuous family of new 
many-body problems, labeled by a flow parameter t. In the limit t ^ cx) it sometimes happens that the resulting 
problem simplifies and can be solved in a controlled manner. In order to obtain a non-trivial RG fixed point for 
t — > 00 which is characterized by a non-zero anomalous dimension 77, it is crucial that the fields are properly rescaled, 
as discussed in detail by Bell and Wilson For example, by applying the exact functional RG to ^'^-theory |p^ , |l6| , 
it is easy to see that with out the rescaling the RG flow does not reproduce the Wilson-Fisher fixed point below four 
dimensions, see Sec.LVC. We therefore believe that the RG calculations of Refs. [0-^ are incomplete, and should 



be augmented by the usual rescaling steps 2 and 3 given above. Note that the rescaling is explicitly included in the 
one loop momentum-shell approach advanced by Shankar [|l^ , so that one may wonder why the exact functional RG 
used in Refs. does not include the rescaling. We shall try to clarify this point in this work by showing how 

to include the rescaling of the degrees of freedom in the exact functional RG for fermionic many-body systems in 
arbitrary dimensions. We then point out several advantages of working with the rescaled version of the exact RG flow 
equations. 
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II. MODE ELIMINATION FOR THE LEGENDRE EFFECTIVE ACTION 



In this section we present a brief derivation of the exact flow equation describing the change in the Legendre effective 
action T{if, cp} due to the eUmination of the degrees of freedom with momenta and frequencies in a suitably defined 
sheU. Keeping in mind that the Legendre effective action is the generating functional of the one-particle irreducible 
correlation functions JTst , an expansion of T{ip,ip} in powers of the fields (p and (p yields the corresponding flow 
equations for the irreducible n-point vertices. For classical field theories, the flow equation for the Legendre effective 
action has been considered some time ago by Nicoll et al. and more recently by Wetterich and by Morris 
. For non-relativistic fermions the exact flow equation for the Legendre effective action has recently been derived 
in a very general form by Honerkamp, Salmhofer and co-authors [^||, who also included the possibility of symmetry 
breaking. Our derivation given below is less general but more explicit; in particular, we shall adopt as much as possible 
the standard notations from the theory of critical phenomena pjl^. 



A. Assumptions and initial conditions 

We consider a translationally invariant non-relativistic system of interacting fermions with a Fermi surface. 
Throughout this work we assume that the fermions are in the normal state, so that the propagator does not have any 
anomalous components. By Fermi surface we mean the true Fermi surface of the interacting system, which is defined 
as the set of all momenta satisfying 

ek^ -/i-I](kj.,iO) , (2.1) 

where ek is the energy dispersion of the non-interacting system, /i is the exact chemical potential and S(k, iujn) is the 
exact self-energy of the interacting many-body system. Here k is the momentum and itUn is the Matsubara frequency. 
Given an arbitrary k, we may define a corresponding momentum ki? on the Fermi surface via the decomposition pc| ] 

k = kF+VFP, (2.2) 
as shown in Fig.|l|. Here irp — vf/WfI is a unit vector in the direction of the local Fermi velocity 

VF = Vk eklk=k, • (2-3) 
Note that vp is defi ned in terms of the gradient of the bare energy dispersion at the renormalized Fermi surface. The 



solution of Eq.(2.1) can be parameterized as 

kf = nfcF(n) , (2.4) 

where fi is a unit vector in the direction of k^, see Fig.|l|. Note that only for a spherical Fermi surface we may 
identify n = vp, so that ki? and vj? are parallel. In this case we know that k^ is not renormalized as we turn on the 
interactions at constant density The true Fermi surface is then given by ekj? — fJ-o (where /io is the chemical 

potential of the non-interacting system at the same density), so that S(kF, iO) = fJ- — fJ-o-, which is independent of k^. 

For the derivation of the RG flow equations it is convenient to represent the generating functionals of the fermionic 
correlation functions in terms of Grassmannian functional integrals p^ . We assume that initially the Grassmann 
fields with large excitation energies ek — ek^ ^-iid large Matsubara frequencies LOn have been integrated out. Here 
K = (k, iuJn) is a composite label, and for simplicity we ignore the spin degree of freedom. The high energy shell 
can be described by an equation of the form ~ , where VIk is a suitable homogeneous function of ek — ^k^ or 
\ujn\- Possible choices are fix = |wn|; = |ek — ekjrl, or fix = ^oj^ + (ek — ekp)^. Our starting point is an effective 
action of the form 

SMi^} = SlM.i^} + S'^'{i},i;} , (2.5) 

where the free part is given by 



J K 



(2.6) 

Here we use the notation 



2 



K 



-y 



f3V ^ J (27r)^y_^2^ 



where /3 is the inverse temperature and V is the volume, and the right-hand side is vahd for /3 — > cxo and V ^ oo. 
The interaction part is 



(2!)^ Jk' Jk' Jk 



Sk[+K'^,K2+Ki 



■2JK1 

xTf^{K[,K'^-K2,K^)i^K[i^K',^K2^K^+... , (2.8) 
where the eUipsis denotes three-body and higher order interactions, and 

^ (27r)^+M(k - k')(5(t^ - tj') • (2.9) 

An vertices are antisymmetric with respect to the permutation of any pair of the incoming particles and any pair 
of the outgoing particles. In particular, the four-point vertex t''^\k[, K2; K2, Ki) is antisymmetric with respect to 
the exchange K[ ^ K2 and Ki ^ K2. The function 67(e) is a smooth cutoff functi on, satisfying 0^(6) « 1 for 



e > 7 and 67(e) « for e < —7, so that lim^^o 67(E) = 6(e). The term 'S^g{K) in Eq.(2^) is the contribution from 



the high energy fields with CIk ^ £n to the irreducible self-energy, and ~Yi(kp, iO) is a counterterm which takes into 



account that in the free action (2_^) we have subtracted ekj? = /i — I](ki?,iO) from the bare energy dispersion. As 
shown below, this subtraction is crucial to obtain the RG flow of the Fermi surface. In a perturbative approach, the 
above subtraction is necessary to obtain a well-behaved perturbation series [|l9|j23| ]. The counterterm — I](ki?,iO) can 
be determined from a self-consistency condition, which can be imposed order by order in perturbation theory [ p4[ . 
Within a RG approach, the counterterm can be determined a posteriori from the condition that the flow equation 
of t he nr omentum- and frequency-independent part of the irreducible two-point vertex has a fixed point p^ , see 



Sec^VB] for a careful discussion. The energy scale ^0 plays a dual role ||13|: For the calculation of the vertices of the 
initial action S^„{ip,tp} the scale ^0 acts as an infrared cutoff, so that for sufficiently large ^0 the vertices appearing 
in S'^qIi/j, V'} can be calculated perturbatively. On the other hand, for the remaining low-energy degrees of freedom 
plays the role of an ultraviolet cutoff. 



B. Exact flow equation describing mode elimination 



Starting from the effective action S(^g{'ip,'ip} with cutoff ^q, we eliminate all fields with momenta and Matsubara 
frequencies in the regime f ^ ^ S,q ■ The free propagator for the fermionic fields in this regime has in the if -basis 
the diagonal elements 



Ck + ekj 



(2.10) 



Because lip lies by construction on the true Fermi surface of the interacting system, the right-hand side of Eq.(2.1C) 
depends implicitly on the interaction. Note, however, that ekj^ is independent of the flow parameter ^. Diagrammati- 
cally, the elimination of the degrees of freedom in the shell ^ ^ fix ^ ^0 corresponds to contracting all terms generated 



by expanding e "^so 



with the propagators G^^^ given in Eq.(pl^). The connected correlation functions of the 
new theory can be formally represented as functional derivatives of the generating functional G^{J, J} defined by |2^] 



(2.11) 

Here (a, Mb) — J^^ Jj^, a^MxK'bK' , and the factor ( = —1 arises from the antisymmetry of the Grassmann fields 
p2| , p6[ . should be considered as a matrix in if-space, with matrix elements given by 
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(2.12) 



where G°,^„(ii') is defined in Eq.( plo|). 

Differentiating botfi sides of Eq.(2.11) with respect to ^, we obtain an exact flow equation for the generating 
functional of the connected correlation functions. For practical calculations the flow equation for the irreducible 
vertices is more convenient [ pj|jl2lJl^ . To obtain the corresponding generating functional, we perform a Legendre 
transformation. 



c^{^,^}^{^,j) + {j,v)-gi{j,j} 



(2.13) 



where the Grassmann sources J and J have to be considered as functionals of the Grassmann fields ip and ip by solving 
the following equations for J = J{<f', and J = Jiip, </?}, 



SGHJ, J} 



5J ' SJ 
The generating functional of the irreducible vertices is then given by 

ra<p, ^} - c^{^, ^} - (^, [-G^.^J- V) 

and satisfies the exact flow equation [p|j5[|27|| 



(2.14) 



(2.15) 



Jk 

+Gl,M^-Gl,t,rGl,^} 



22;KK 



with initial condition 



(2.16) 



(2.17) 



Here {. . ■}22-kk denotes the lower diagonal element of the corresponding 2 x 2-matrix in X-space, and the functional 
is given by 



c 



C-, 



KK' 



KK' 



(2.18) 



where is defined as the field-independent part of the second functional derivative of F^, so that ZY^{0, 0} = 0. The 
interacting cutoff-regularized propagator is related to the non-interacting one via the Dyson equation, 



Gi..io - [G",4„ - S^] ^ , [G"^^Jjj - (Jy-G^^^^ 



(2.19) 



C. Flow equation for sharp cutoff 



For simplicity we shall from now on work with a sharp cutoff function, 6(e) — lim^_,o 0^(e). Due to translational 
invariance the if-dependence of all vertices is constrained by energy-momentum conservation, so that the matrix 
elements of are 



[^i\ij-KK' = 5vj5KK'[^i{K) - i;(k_F,zO)] , 



(2.20) 



where the term E(ki?,iO) is due to the subtraction in Eq.( p.8[ ). The exact fiow equation (2.16) can then be reduced 
to the following form [p^,p^ 
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We now expand 



c 



^ ta^„ — ek + /i — (if) I. " )22;KK 

iun - ek + A* - (K) 



-CPV / 6{nK - In 
Jk 



Ck + ekf 



(2.21) 



n=Q 
X Sk' 



K'K„- 



Ki 



Ki 



X (fK{ ■ ■ ■ <^K; fKr, ■ ■ ■ fKi , 



(2.22) 



and identify the terms with the same powers of the fields on both sides of Eq. ( ^.21 ) . Note that in the normal state only 
the even vertices are non-zero. In this way we obtain the RG flow equations for the unrescaled irreducible 27i-point 
vertices F^^"'' {K'l , . . . , K'^; Kn, ■ ■ ■ , Ki). We now explicitly give the exact flow equations for the vertices r^"-* , r^^-* , 

and T^^\ For a two-loop calculation one needs the flow equation for the six-point vertex T^^\ which is given in the 
Appendix. 



1. Free energy 



The interaction correction to the free energy is obtained from the last term in Eq.(2.21) 

d^vf^ - / 5{nK - C) In 

JK 



iuj„ - £k + /-^ - ^d^<) 



(2.23) 



2. Self-energy 



Comparing the terms quadratic in the fields on both sides of Eq. ( |2.2l| ) and using the fact that by construction the 
two-point vertex is related to the subtracted irreducible self-energy via 



Tf{K-K) = -[^^{K) - S(ki.,zO)] 



we obtain 



(2.24) 



In, iujn' - Ek' + A* - ^iiK') 

A graphical representation of this equation is shown in Fig||. 

3. Four-point vertex 
The flow equation for the irreducible four-point vertex is 



(2.25) 
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K 



J71 - Ck + - '^d^) i^n' - fk' + /-i - ^d^') 



1 

2 

+c 

+ 



r'f\K'„K'; K, K,)Tf\K[,K- K', K^) 



(4)/ 



K'=Ki+K2-K 

I 

. K'=K+Ki-K[ 
K'=K+Ki-K' 



(2.26) 



This equation is shown graphically in Figj^. The first term in the curly braces is contribution from the BCS-channel, 
while the last two terms are the zero-sound contributions, usually abbreviated by ZS (second term) and ZS' (third 
term) Note that both zero-sound terms have to be retained in order to preserve the antisymmetry of the 

four-point vertex psl . The exact flow equation for the six-point vertex is rather lengthy and is given in the Appendix. 



III. EXACT FLOW EQUATIONS DESCRIBING MODE ELIMINATION AND RESCALING 

So far we have derived exact flow equations for the irreducible vertices describing the elimination of the degrees 
of freedom. Within the one-loop approximation, it is sufficient to set r^^"''' = for n > 3, and to ignore interaction 
corrections to the propagators in internal loops. The resulting truncated flow equation for the four-point vertex in two 
dimensions has been analyzed numerically by Honerkamp et al. |^,^. A similar numerical analysis of an equivalent 
truncated flow equation has been performed by Halboth and Metzner Both groups found that the one-loop 

flow of the four-point vertex eventually diverges at a finite scale, where the perturbative RG breaks down. Physically 
the runaway flow has been interpreted in terms of some incipient instability of the normal metallic phase. Here we 
would like to point out that in principle there is another interpretation of this runaway flow, namely the existence 
of a scale-invariant non- Fermi liquid fixed point, which is characterized by a finite anomalous dimension. Below we 
argue that the mode elimination RG transformations of Refs. fll-^ cannot detect such a fixed point, because these 
equations do not take into account that the degrees of freedom should be properly rescaled to reach a fixed point 
with a finite anomalous dimension Assuming the existence of such a fixed point, only the rescaled version of the 
exact flow equations given below would detect it, while the pure mode elimination RG used in Refs. would still 
exhibit a runaway flow to strong coupling. 

Because the following rescaling procedure depends crucially on the existence of a Fermi surface, we now explicitly 
set the statistics factor ( = —1. We also take the limits of infinite volume {V 00) and zero temperature (/3 00), 
so that momenta and frequencies become continuous variables. 



A. Scaling variables 



Instead of the momentum k and the Matsubara frequency ilu, we now label the degrees of freedom by the direction 
n of ki? (see Fig.|l|) and by the dimensionless variables 

VFP VF-(k-kF) 

<?= — = ' '^"i' ^ ' 

so that 

]i. = 'kp+-vp — q = nkp(h) + vp — q. (3.2) 

Vp Vp 

Furthermore, instead of the flow parameter ^, we introduce the dimensionless logarithmic flow parameter 

t = -ln(e/eo), (3.3) 

and define the dimensionless dispersion 
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where 



(3.4) 



(3.5) 



is an irrelevant coupling which meas ures the leading deviation from linearity in the energy dispersion in the direction 
normal to the Fermi surface. In Eq.(3.4) we have expanded the energy dispersion around k — k^?, 



Ck — ekf 



• (k — kj; 



(k-kf^)- 
2 m 



(3.6) 



Note that in general "vp and m depend on n. Using Q — (n, q, ie) instead oi K — (k, iuj) as integration variables, we 
have for /3 — > oo and V oo 



K 



dSn 

Sd 



dqJ{h,q) 



2^ 



(3.7) 



where dS^ is a surface element and So is the surface area of the unit sphere in D dimensions, and J(n, q) is a 
dimensionless Jacobian associated with the transformation k (n, g). For convenience we have pulled out a factor 
of z^oC^ in Eq-(^-7), where i/q is the density of states at the Fermi surface, 



1^0 



dk 



With this normalization J(n; q) is dimensionless. In particular, for a spherical Fermi surface 



J{h,q) = (l + ctq) 



D-l 



(3.8) 



(3.9) 



B. The scaling form of the irreducible vertices 



The proper definition of the dimensionless scaling form of the irreducible vertices follows partially from dimensional 
analysis, and partially from aesthetic consi derat ions (such as the requirement that numerical prefactors should be as 
simple as possible). Given the expansion (2.22) of the generating functional of the irreducible vertices in powers of 
the dimensionful fields ipK, we substitute 



zi_ 



1/2 



(3.10) 



where the wave-function renormalization factor is related to the irreducible self-energy E^(k, iw) at scale ^ as 
usual, 



1 



9S;(kF, 



(3.11) 



The dimensionless fields (pQ should be considered as functions of the scaling variables Q. 



1. Free energy 

Due to the Jacobian associated with the rescaling of fermionic fields in the functional integral, the field-independent 
part r^°^ of the generating functional T^{<p, (p} picks up an additive term, so that after rescaling the correction to the 
free energy is 

f (0) ^ p(o) ^ /■ e(^_f]^)inZf . (3.12) 
Jk 

In classical statistical mechanics a contribution analogous to the second term has been discussed by Wegner and 
Houghton 
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2. Two-point vertex 



We write the exact propagator in the following scaling form, 

„ . . Z^-fkp VF ■ (k - kp) iuj 



(3.13) 



where fi should be considered as a function of k, as given in Eqs.( |2.2| ) and (2^). Note that our exact RG equations 
describe the flow of the dimensionless scaling function 



GtiQ) = Gtin,q,ie) = J^G^.^^ (nfcF(n) + vf— g, i^e) ■ 
Introducing the dimensionless scaling form of the (subtracted) irreducible two-point vertex, 



f['\Q) = ^Tf\K) = [J^iiK) S(kF,zO)] 



and the scaling form of the inverse propagator 



we obtain 



rt{Q) = [iw - Ek + M - ^^(k, it^)] 
= Zf [z6-ef(9)]+ff^(g), 

Q(SIq - 1) - Q{nQ - e') 



GtiQ) 



niQ) 

e(e* > fig > 1) 



where 



and 



„ , , I 1 if a;2 > a; > xi 



n 



K 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



Note that for the choice VIk = |ek — I '^^ obtain (1q = \£,i^{q)\ ~ \q\ to leading order. The initial condition at t = 
implies 



ro(Q) = T^N-ek + M"S^oW] 
Furthermore, without interactions Zt ~ 1 and '{Q) — 0, so that rt{Q) = ie — ^"(g) 



(3.20) 



3. Higher order vertices (n > 2) 



The dimensionless scaling form of the higher order irreducible vertices follows from our definition (3.10) of the 
dimensionless fields. 



-^t HQij ■ ■ ■ I Qn'j Qni ■ ■ ■ 1 Ql) 



,n— 1 (-n — 2 



Ff"\x;,...,<;if„,...,ifi) 



1/2 



(3.21) 
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C. Flow equations for the rescaled vertices 

1. Free energy 

Defining the interaction correction to the free energy per Fourier component, 



r(o) 



and using Eqs.(2.23) and ( 3.12| ), we obtain the exact fiow equation for ft, 



dtft ^ ft^ J '^'^t 
Here ry" is the flowing anomalous dimension, 



Ml 



(3.22) 



(3.23) 



(3.24) 



Eq.( |3.23 ) is the fermionic analog of the corresponding flow equation for the free energy in (/)^-theory, see Eq.(4.6) of 
Ref.W. 



2. Two-point vertex 



From Eqs.(2.25) and (3.15) we find the flow equation for the dimensionless subtracted two-point vertex, 

dttf\Q) = {l-i^f-Q-dQ)Tf\Q) 



where we have introduced the notation 



Gt{Q')t['\Q,Q';Q',Q) 



Q ■ Oq ^ qdq + ede , 
Sing - 1) 



Gt{Q) 



niQ) 



(3.25) 

(3.26) 
(3.27) 



3. Four-point vertex 

The flow equation for the rescaled irreducible four-point vertex follows from Eqs.( 2.26| ) and ( 3.21 ) 
dtf[''\Q[,Q',;Q2,Qi) = 



E 



f^'^(Q'i,Q2;Q2,Qi) 



Gt{Q)f['\Q[,Q'^,Q;Q,Q2,Qi) 
Gt{Q)Gt{Q') + Gt{Q)Gt{Q' 
f^t\Q'i.Q'2,Q',Q)T^t\Q,Q':Q2,Qi) 

ff\Q',,Q'-Q,Q,)ff\Q'^,Q-Q',Q2) 
'f^\Q'2, Q'; Q, Qi)ff\Q{, Q; Q', Q2] 



K'=Ki+K2-K 
K'=K+Ki-K[ 

K'=K+Ki-K' 



(3.28) 
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Here K', Ki and K[ should be considered as functions of the dimensionless scaling variables introduced in Sec. HI A , 
for example K = [nkp + '^f:^Q, Due to the non-linearity of this transformation, the explicit expression of Q' in 
terms of Q and Qi — {hi, qi, e^) is rather complicated. For example, let us calculate Q' — {n' , q', e') in the zero-sound 
contribution involving K' = K + Ki — K[. The energy component is simple, e' = e + ei — e\, but for n' and q' we 
obtain 



n + ni - -t- ct(ng -f fiigi - h'^q[) 
|n -f ni - + Ct{nq + ftigi - n[q[)\ 



(3.29) 



q' = — [|n + ni - n'^ + ct{hq + hiqi - h[q[)\ - 1] 



(3.30) 



where c* = Coe */{mvp), see Eq.(3.5). For simplicity we have assumed a spherical Fermi surface, so that k,; 
hikp{l + ctqi). Note that the corresponding rescaled energy dispersion is 



{q')=q + h-{hiqi-h[q[] 



fi ■ (fii - h[) 

Ct 

2ct 



+ {hi - h[) ■ {hq + hiqi - h[q[) . 
The rescaled flow equation for the six-point vertex is given in the Appendix. 



(3.31) 



IV. ADVANTAGES OF RESCALED FLOW EQUATIONS 



In this section we discuss several properties of the rescaled flow equations of Sec. Ill, and argue that for practical 
calculations it may be advantageous to use the rescaled flow equations instead of the unrescaled equations discussed 
in Sec.||. 



A. Screening and flow of the density-density correlation function within RPA 



As pointed out in Refs. [||-D, the unrescaled flow equations predict pathological RG flows for physical quantities 
which are determined by degrees of freedom in the immediate vicinity of the Fermi surface, such as uniform suscep- 
tibilities: for any finite infrared cutoff ^ uniform susceptibilities are not renormalized at all, while at ^ = their RG 
flow exhibits a discontinuity. Although for a Fermi liquid one may use the Fermi liquid relations between the uniform 
susceptibilities and the quasi-particle interactions to obtain flow equations for the susceptibilities at flnite cutoff ^ 
1^,^, it would be better to calculate uniform susceptibilities entirely within the framework of the exact RG, without 
relying on the assumption that the system is a Fermi liquid. 

To illustrate the above point, let us consider the screening of the effective interaction, which is closely related to 
the density-density correlation function. The screening problem has also been considered by Shankar (see Appendix 
A of Ref. 1^^) within the flcld theory version of the RG, and by Dupuis within the conventional Kadanoff- Wilson 
RG scheme for Fermi liquids. 



Consider the flow equation (2.26) for the irreducible four-point vertex. Suppose we start from a bare interaction at 
scale ^0 of the form 



rfJ{K[,K!,; K2,Ki)^ U,{Ki ~ K[) - U,{Ki - K'^ 



(4.1) 



Let us now calculate the flow of the effective interaction within the random-phase approximation ( RPA ) , where the 
BCS-contribution is ignored and only those zero-sound terms are retained which preserve the form (4.1). Setting for 
simplicity Ki = K[ and K2 — K2 (this is the combination appearing in the flow equation ( 2.25| ) for the self-energy) 
we approximate the irreducible four-point vertex at scale ^ < ^0 by 



Tf{Ki,K2;K2,Ki) « -/^(i^i -ifa) , 



(4.2) 



where we have assumed /j(0) = 0. In this approximation we obtain from Eq.(2.26) 
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where the change of the polarization is given by 

' Si^K - OG^^UK + Ki - K2) 



(4.3) 



K 



1^71+711-712 - ek+ki-ka + ^ — Sj(X + Ki - K2 



(4.4) 



Within the RPA we ignore the self-energy corrections to the propagators in Eq.(4.4). Using the cutoff function 
^^K ~ |ek — Ck^lj we obtain for |ki — k2| <^ kp 



tl^{Ki - K2) ^ 26(0 



dk 



(2^)U'5(ek-6kJ 
vp ■ (ki - k2) 



l(w„i - UJn^ ) - V_F • (ki - k2) 



(4.5) 



Obviously, ri^(P) = for any ^ 7^ 0, so that the RPA interaction f^{P) is not renormalized for any finite ^. Due to 
the factor 6{^) in Eq.(4^) the RG flow of the RPA interaction f^{P) (and hence the flow of the polarization and the 
cornpressibility) is discontinuous at ^ = 0. This discontinuity is smoothed out if one works at a finite temperature 
p9| , p0| ; but there is a way to avoid discontinuous flow equations at zero temperature: the RG flow equations generated 
by our rescaled version of the exact RG remain continuous even at T = 0. 

We now explicitly show this for the RG-flow of the rescaled density-density correlation function. Within the RPA 
this the rescaled four-point vertex is approximated by 



Tt\Qi,Q2;Q2,Qi) ~ 'FtiQ,,Q2) , 

where the dimensionless function Ft{Qi, Q2) satisfies the flow equation 

[dt - Qi ■ dQ, - Q2 ■ dQ,]F-\Q,, Q2) = -n,(Qi, Q2) , 



which is the rescaled analog of Eq.(4.3). Here 



nt(Qi,Q2) 



GtiQ)GtiQ') + Gt{Q)Gt{Q') 



K'=K+Ki-K2 



(4.6) 



(4.7) 



(4.8) 



where Q' = Q'{Q,Qi,Q2) is defined as a fu nction of Q = (n,q,ie), Qi — {hi, qi,iei) and Q2 — {ii2,q2,'i£2) via 
K' = K + Ki — K2, as explained in Sec. Ill C. Ignoring again interaction corrections to the propagators and working 
with the momentum shell Q,q = |ek — ek^l/C — IC"!?)!' have 



Gt{Q) 



Q{e* > lef (g)| > 1) 

*e-^f(g) 



(4.9) 



GtiQ) 



Substituting these expressions into Eq.(4.8) and performing the integrations over q and e we obtain 

dSn 6(1 + e* > n • qi2 > 2)n • qi2 



nt(Qi,Q2) = -4 



Sd (n • qi2)2 + (ei - 62)^ 



where we have used Eq.(3.31) to simplify ^" (g') for small |ki — k2 



(<?') ~ <? + n ■ qi2 , qi2 = nigi - n2(72 



ni - n2 

ct 



(4.10) 



(4.11) 



(4.12) 
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The crucial point is now that, unlike Eq.(4.5), the right-hand side of Eq.(4.11) is a non-singular function of the flow- 
parameter t. In fact, in dimensions D > 1 the 0-function discontinuity of the integrand disappears after the angular 

integration, so that Ilt{Qi, Q2) is a continuous function of t. But even in 13 = 1, where the angular integration should 
be replace d b y a summation over the two discrete Fermi points, there is no discontinuity in the solution of the flow 
equation (Ej). 



It is instructive to elaborate a little bit more on the case Z3 = 1, where we obtain from Eq.( 4.11 ) 



nt(Qi,Q2) 



2<^ni,n2 



|gi -g2|e(l + e* > |gi -g2| > 2) 
(gi - 92)2 + (ei - £2)' 



(4.13) 



In deriving this expression we have assumed Ct\qi\ ^ 1, and have neglected the terms with fii = — n2, which vanish 
exponentially for t — > 00. Note that the right-hand side of Eq.(4.13) has discontinuities, but no (5-function singularity. 



in contrast to Eq.(4.5). It follows that the flow of the physical density-density correlation function is continuous. To 
see this, let us explicitly solve Eq.(4.7). Using the method described in Ref. |jl^, the solution with the correct initial 
condition is easily obtained. 



e '^zei; 112,6 



92, e 



Substituting Eq. ( [4. 131) into Eq. ( [4.14| ) we finally obtain in D = 1 



Jlt{qi - 92, iei - i£2) 



(4.14) 



(4.15) 



where 



n((g,ie) 



^^[(|g|- 2)6(1 + e*>M>2) 
(2e*-|g|)e(2e*>M>l + e*)] . 



(4.16) 



Note that the function xQ{x) is continuous, so that the right-hand side of Eq.(4.16) is indeed a continuous function 
of t. The physical density-density correlation function is then 



rr ^ ■ \ ^ f VFP 

-\p\ivF\p\~20ei^ + ^o>VF\p\>20 



{vfpY + 

\pm^-VF\p\)Q{2io>VF\p\>i + ^o) 



{vfpY 



- W 



Taking the limits ^0 00 and ^ 



we obtain the well-known result 

„2 



no(p,«w) 



VF 



P 



TT {vfpY + 



(4.17) 



(4.18) 



From the first line in Eq.( 4.17 ) it is clear why our rescaled flow equations are more suitable for the calculation of 
the susceptibilities than the corresponding unrescaled flow equations: Because the dependence of the density-density 

correlation function on the infrared cutoff appears in the scaling function lit via the ratios p/^ and w/^, 

the result for uniform susceptibilities depends on the order in which the limits p — > and ^ — > are taken. Our 

rescaled RG equation directly yields the scaling function tit (^-^j where this problem does not arise. Let us give 

an alternative explanation for this difference: in the unrescaled flow equations the momenta and frequencies are held 
constant as the infrared cutoff is reduced. Any flxed distance from the Fermi surface is therefore magnified on the 
scale of the reduced infrared cutoff, so that momenta that are initially close to the Fermi surface are mapped onto 
new momenta which, on the reduced scale, appear further away from the Fermi surface. On the other hand, in the 
rescaled fiow equations the momenta and frequencies are scaled down together with the infrared cutoff, so that the 
degrees of freedom that are initially in the vicinity of the Fermi surface are mapped onto coarse grained degrees of 
freedom that remain close to the Fermi surface. 
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B. The Fermi surface as a RG fixed point manifold 



As already emphasized in Sec. II A , to obtain well-behaved scaling properties we should expand the bare energy 
dispersion around the true Fermi surface of the interacting many-body problem, which a priori is not known. We 
now show how the shape of the Fermi surface can be calculated a posteriori from the fixed point equation for the 
momentum- and frequency-independent part of the irreducible two-point vertex, 



rr(n,o,iO). 



(4.19) 



Note that this is a relevant coupling function (labeled by the direction fi) with scaling dimension -1-1. There are 
two marginal coupling functions associated with the two-point vertex. One of them is the usual wave-function 
renormalization given in Eq.(3.11). In terms of our rescaled two-point vertex it can also be written as 



d{ie) 



(4.20) 



e=0 



The other marginal coupling can be taken to be the dimensionless renormalization factor of the Fermi velocity 



9ff)(n,g,zO) 



dq 



(4.21) 



g=0 



The RG flow of the shape of the Fermi surface of strongly correlated electrons in reduced dimensions have recently 
been considered by several authors ]l7|,p|,^ 31 1. However, the interpretation of the renormalized Fermi surface in 
terms of a RG fixed point manifold has only been emphasized by Shankar p^ ], who used the conventional one- loop 
momentum shell technique, including the rescaling step. In contrast, the authors of Ref. ||,^,|l|] work with unrescaled 
flow equations, in which case the coupling corresponding to Eq.( 4.19| ) does not flow to a RG fixed point, so that the 
interpretation of the Fermi surface as a RG fixed point manifold is obscured. 

To study the shape of the Fermi surface within the exact RG, let us retain only relevant and marginal couplings in 
the two-point function. In this approximation 



(4.22) 



It is convenient to define the Fermi surface momenta \<iF,t of the system with cutoff ^ = ^oe * in analogy with Eq.( ^.l[ ), 

ek^^ =^-S^(kF^t,iO) , (4.23) 



so that the true renormalized Fermi surface is k^^ = limi^ookF.t, see Eq.(2.1). In order to calculate k^^^t we use 
the fact th at by construction q in Eq.(4.22) is measured relative to the renormalized k^^, so that we obtain from the 
definition (4.23), up to irrelevant terms. 



kp t = k/ 



•vp- 



(4.24) 



Setting t — Q we obtain for the difference between the fully renormalized Fermi momentum k^? and the Fermi 
momentum \i.F,o of the model with cutoff 



(4.25) 



It should be ke pt in mind t hat v j? is defined in terms of the gradient of the hare energy dispersion at the renormalized 
\^F■ From Eqs.( [3.15 ) and (4.19) we find 



CoAj? = ^o"P(kF,iO) - %(kf ,iO)] , 



(4.26) 



which allows us to reconstruct the counterterm I](ki?,zO) as follows: Suppose that we have adjusted the initial 
condition /i" such that for i — s- cx) the relevant couplings /2" approach a fixed point. Due to the relevance of /i", 
we expect that this requires some fine tuning of the initial /i" , so that on the critical manifold (i.e. the manifold in 
parameter spac e that flows into a RG flxed point) /i" becomes a function of the other relevant and marginal couplings. 
From Eq.( fl.26 ) we then obtain the counterterm E(ki?, iO) as a function of these couplings. Note that if we adjust the 
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initial conditions such that fif remains finite for t — > oo, then Eq.(4.24) guarantees that ioi t ^ oo the flowing k.p,t 
indeed approaches the true Fermi surface. Thus, the problem of constructing the renormalized Fermi surface can be 
reduced to the problem of finding the RG fixed point o f the flow equation for the relevant coupling function jlf. 

The above discussion relies on the approximation ( |4.22| ) for the inverse propagator, which is only justified for a 
Fermi liquid. In principle it is also possible that the system flows to a non-Fermi liquid fixed point, characterized by 
some new renormalized Fermi surface. As pointed out by Anderson ||3^ , in this case there might be some subtleties 
related to the subtraction of the counterterm S(ki?, iO). 

It is also worth emphasizing the following point: In order to reach a RG fixed point, it is crucial that the bare 



energy dispersion is expanded around the fully renormalized Fermi surface ki?, defined in Eq.(2.1). If we had chosen 
the bare Fermi surface \<iF,o as the reference for the expansion, then the RG flow of the relevant coupling /i" would 
exhibit a runaway flow, i.e. — > o o for t oo. To see this, suppose that we expand around the initial Fermi 



surface at scale ^o- Then the relation (4.24) between the renormalized and the bare Fermi momentum is replaced by 



kf,«^kf,o + Vf.o ^°^ If , (4.27) 

where vp^Q = VkCklkpg- Given that lip.t approaches for large t a value that is different from the bare kj^^Oj we see 
from Eq.(4.27) that fif must necessarily diverge as e* (assuming that vf remains finite). This runaway flow indicates 



that we have expanded around the wrong Fermi surface. This phenomenon is well known from the usual theory of 
critical phenomena, where the fixed point manifold is discrete: In this case a runaway flow indicates the existence 
of some new fixed point. What is new here is that the fixed point manifold is a continuum, which can therefore 
continuously change under RG transformations. To avoid a runaway flow, we have to expand all quantities from the 
beginning around the true Fermi surface, which can be determined a posteriori from the requirement that the RG 
flow indeed approaches a fixed point. 



C. Rescaling and the possibility of non- Fermi liquid fixed points 

Long time ago Bell and Wilson ^ pointed out that the field rescaling (see step 3 in SecJ|) is necessary in order to 
obtain a non-trivial RG fixed point, characterized by a non-zero anomalous dimension 77. In fact, the RG fixed point 
conditions can be viewed as some sort of non-linear eigenvalue equations, which contain 77 as an adjustable parameter 
p3| . From this point of view it is not surprising |l5| that solutions of the fixed point equations exist only for certain 
values of 77, depending on the interaction, the dimensionality, and the symmetries of the system. We believe that the 
above statement is also true for fermionic many-body systems, so that the unrescaled flow equations used in Refs. 
|^|-^ cannot be used to detect non-Fermi liquid fixed points. 

Because within a one-loop approximation the anomalous dimension vanishes, at the first sight it seems that at this 
level of approximation it is sufficient to work with a pure mode-elimination RG This is not necessarily true, 

however, because even at the one-loop level the rescaling of momenta and frequencies (see step 2 in Sec.^ can be 
essential to obtain a non-trivial RG fixed point. For example, consider (/)^-theory slightly below = 4 dimensions. In 
this case there are two relevant couplings, namely the momentum-independent part of the irreducible two-point vertex 
fit, and the momentum- independent part of the irreducible four-point vertex gt- The one- loop RG flow equations for 
suitably defined ||l^ dimcnsionless couplings are well known ||l^ 

2(1 + m) 

dfgt = (4 - D)g, - . (4.29) 

2(1 -I- ^itY 



For < 4 — D ^ 1 these equations have a non-trivial fixed point at 



li^^D) , ^*«-i(4-i?), (4.30) 
which arises from a bala nce b etwee n the competing effects due to the scaling terms (i.e. the first terms on the right- 



hand si des o f Eqs.( 4.28 ) and (4.2£)) and the mode elimination terms (i.e. the second terms on the right-hand sides 
of Eqs.( 4.28 ) and ( 4.29D ). Obviously, a RG transformation which includes only the mode elimination would miss the 



Wilson-Fisher fixed point. 
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Finally, we point out that the rescaling of momenta and frequencies plays also an important role to accelerate the 
flow of the irrelevant couplings to negligibly small values. For example, the constant part of the rescaled irreducible 
six-point vertex vi^^Q'-, , Q9, Q'ri; Q3, Q 21 Qi) is irrelevant with scaling dimension —1. This —1 appears as the first 



term on the right-hand side of Eq.(A.4), and causes the size of the rescaled six-point vertex to diminish exponentially 



as we iterate the RG (as long as t his d ecrease is not overwhelmed by some anomalous dimension). On the other hand. 



the corresponding flow equation ( |A.lD for the unrescaled six-point vertex does not exhibit a similar suppression due to 
a negative scaling dimension, so that there is no guarantee that the unrescaled six-point vertex really becomes small 
as we iterate the RG transformation. 

V. SUMMARY AND CONCLUSIONS 

In this work we have constructed exact functional RG equations for non-relativistic fermions, taking not only the 
mode elimination into account, but also the rescaling of momenta, frequencies, and fields. The exact flow equations 
given in this work are valid in arbitrary dimensions and describe the RG flow of any translationally invariant fermionic 
many-body system in the normal state. We emphasize that these equations are exact and describe the RG flow of the 
irreducible correlation functions of the many-body system for all momenta and frequencies. Of course, it is possible 
to extract the usual perturbative RG /^-function describing the flow of the momentum-independent part of the four- 
point vertex from these flow equations. To obtain the /3-function at the two-loop order, it is sufficient to truncate the 



inflnite hierarchy of flow equations by setting F^^"^ ~ for n > 4. In Sections 11 C, [II C and in the Appendix we 
have explicitly written down all expressions necessary for performing a two-loop calculation for any £)-dimensional 
translationally invariant normal fermionic many-body system. These equations are also valid in Z) = 1, where two- 
loop calculations are usually performed using the field theory- method 1^^. Note, however, that in this case there are 
some subtleties which are not fully understood [ ^Sf . 

We have formulated the exact functional RG for fermionic systems in close analogy with the functional RG for field 
theories and statistical mechanics problems. We have tried to emphasize this analogy by adopting as much as possible 
the notation used in statistical mechanics I^Jl^. A crucial step is the introduction of dimensionless scaling variables 



and rescaled vertices, which describe the scaling towards the true Fermi surface of t he sy stem. Finally, in Sec. IV we 



have shown that the rescaled flow equations for the irreducible vertices given in Sec. [II C have several advantages as 



compared with the unrescaled flow equations discussed in Sec.lIC 



1. Flow of susceptibilities. The unrescaled exact RG equations used in Refs. yield discontinuities in flow 
equations for quantities which are dominated by degrees of freedom in the immediate vicinity of the Fermi 
surface, such as the density-density correlation function at long wavelengths. In contrast, our formulation of the 
exact RG including rescaling yields flow equations which are continuous functions of the flow parameter. 

2. Renormalization of the shape of the Fermi surface. The calculation of the shape of the renormalized Fermi 
surface can be reduced to the solution of the fixed point equation for the properly rescaled irreducible two-point 
vertex at vanishing external momentum and frequency. Without rescaling, the corresponding flow equation does 
not have a fixed point, so that the associated relevant coupling exhibits a runaway flow. 

3. Non-Fermi liquid fixed points. If the system has a scale invariant non-Fermi liquid fixed point characterized by 
a non-zero anomalous dimension ry, then the flow generated by RG transformations which omit the rescaling 
steps 2 and 3 listed in Secjj cannot detect this fixed point. Instead, the pure mode elimination RG calculations 
of the type used in Refs. would generate a runaway-flow to strong coupling in this case. 

We suspect that our rescaled flow equations will also turn out to be advantageous for numerical calculations, because 
they describe the RG flow of the relevant scaling functions. The numerical analysis of these equations is a rather 
tedious task which is beyond the scope of this work. 

Our flow equations also offer a new approach to study interacting fermions in D = 1, where the normal metallic 
state is known to be a non-Fermi liquid. For these systems the RG /3-function is often calculated within the field 
theory formulation of the RG advanced by Solyom ||3J]. However, till now there exists no RG calculation of the 
full momentum- and frequency-dependent single-particle Green function of interacting fermions in Z? — 1. From 
bosonization it is known that in momentum-frequency space the single-particle Green function exhibits interesting 
features such as power-law singularities, which are a manifestation of anomalous scaling and spin-charge-separation 
. We are currently investigating whether the full spectral function of interac ting fermions in D = 1 can be calculated 



by means of an approximate solution of the rescaled functional RG equation (3.25) for the two-point vertex. We have 
preliminary evidence [ p7| that at weak coupling this can indeed be done, and that for the Tomonaga-Luttinger model 
the resulting spectral function compares quite well with the exact result obtained via bosonization. 
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Finally, we would like to address the question whether possible scale-invariant non-Fermi liquid fixed points in 
dimensions D > 1 are accessible within the exact RG. Certainly, detecting such a fixed point requires at least a 
two-loop calculation, because at the one-loop level the wave-function renormalization remains unity. Considering 
the complexity of numerical analysis of the one-loop equations describing the RG flow without the rescaling steps 
[^-|| , the direct numerical analysis of the two-loop RG flow equations including rescaling seems to be a rather difficult 
task. Suc h a calculation would require the numerical analysis of the flow equation for the six-point vertex given in 
Eq. (A.4). Note, however, that in the special case of a square Fermi surface Binz, Baeriswyl and Dougot |Q have 
recently presented an analytic study of the one-loop RG flow in the vicinity of the dominant instabilities. Let us 
emphasize again that — > at a non- Fermi liquid fixed point, so that the inclusion of the flow of the wave- function 
renormalization is crucial to reach such a fixed point. 

For Fermi surfaces with a special geometry it is possible to analyze the RG flow analytically [ ^J39| . In particular, 
very recently Ferraz 1 39 1 performed a field-theoretic two- loop RG calculation for a two-dimensional Fermi system with 
a truncated Fermi surface, consisting of flat and curved pieces. He succeeded to calculate the entire single-particle 
spectral function. Interestingly, he found non-Fermi liquid behavior for all points on the Fermi surface, including 
the curved pieces. The calculations presented in Ref. [ p9| are quite encouraging and support our point of view that 
higher-dimensional non- Fermi liquid fixed points are accessible within a two-loop truncation of the Wilsonian RG 
equations presented in this work. 
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APPENDIX: FLOW EQUATIONS FOR THE IRREDUCIBLE SIX-POINT VERTEX 



For a two-loop calculation we need the fiow equation for the irreducible six-point vertex |jT^,|l^,^ . The unrescaled 
flow equation reads (see Fig. ^ and recall that for fermions C = —1) 



-c 

+3 



d^rf\K[,K':„K',;K3,K2,Ki) = 

^^^^ "^Wy^«'' (KlK'^KiK; K, Ks, K,,K^) 



^ icOn — Ek + M ~ (K) 



X <^ An 



(l',2'),3' 



(3,2) 



r'"' {K[,K'2-K',K)vf\K,K\K'^;Ki,K2,K^] 



vf\K[,K'2,K',-K',K, K,)vf {K,K' -Ki.K^) 



K' = K2 + K-i-K 



+9C 



+9 



iiOn - Ek + M ^ '^((K) i^n' - Ek' + A* ~ ^{(-K'') 



X-4(l/_2'),3'-4(2,l),3 



rf ' (K[,K'2, K'; K, K2 , K,)V\'^ {K'-, ,K;K',Kz) 



K'=K+K' -K3 



S{nK^OGi,UK')G^,^,{K") , G5,5o(^)5(t^A-' -OG«„(i^") 



iuj„ — Ek + M ^ S«(-?^) 
G^,^,{K)G^,^,{K')mK" - 



iuj„" - Ck" + M - ^diK" 



iLj„' - Ek' + M - S^(if') 

-4(l',2'),3'-^3,(2,l) 



rf\K[,K'2;K,K')r';>{K'„K'-K", K:i)rY'' {K" ,K'-K2,K, 



K"=K2+Ki-K 
K'=Ki+K2+K3-K' -K 



o.. I 5{nK-0G^,^„{K')G^,^„{K") 

ILJn - Ek + M - ^^(A) 

rf {K[,K'- K, K,)rf^ {K'2, K"- K', K2)Tf {K'„ K; K",Ks) 



K" = K'+K~K3 
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(A.l) 



Here the anti-symmetrization operators ^1.2,3 and ^1,(2.3) are defined as follows, 

^1,2,3/(1, 2, 3) = i [/(I, 2, 3) + /(2, 3, 1) + /(3, 1, 2) 

-/(3,2,1)-/(2,1,3)-/(1,3,2)] , 
A,(2,3)/(l, 2, 3) = ^2,3)4/(1. 2, 3) = i [/(I, 2, 3) - /(2, 1, 3) - /(3, 2, 1)] 



(A.2) 
(A.3) 



Given a function /(I, 2, 3) that is already antisymmetric with respect to the pair (2, 3), the function ^i,(2,3)/(l, 2, 3) is 
a totally antisymmetric function. Note that the combinatorial factors in front of the terms involving the combinations 
r(4)r(6) and r(4)r(4)r(4) are precisely the same as in the corresponding flow equation of (/)^-theory, see Eq.(4.20) of 
Ref. [0. 

The rescaled version of the above flow equation is (setting now C = — 1) 



ftrr (Qi, 02,03 ;Q3,02,Qi) 



-1 



E 



+ Q\ ■ do' + Qi 



ff'(0'i, 02,03; 03,02,00 



Gt(Q)ff (O'l, O2, 03, 0; 0, 03, 02, Oi) 



Gt(0)Gt(0') M(l'.2').3 



+-4i,(3,2) 



f f '(O'l, O2; 0', 0)rr(0, 0', O3; O3, 02, Qi] 



=,(6) 



f r(0'i, O2, 03; 0', 0, Oi)f f '(0, 0'; O3, 02) 



+9 / [Gt(0)Gt(0') + Gt(0)Gt(0')] 



X-4(i'^2'),3''4(2,l),3 



f r(0'i, O2, 0'; 0, 02, Oi)f f '(O3, 0; 0', 



-9 / [Gt(0)Gt(0')Gt(0") + Gt(0)Gt(0')Gt(0") + Gt(0)Gt(0')Gt(0")] Ai',2').3' -43,(2.1) 

f r'(0'i, O2; 0, 0')rr'(03, 0'; 0", 03)rf '(0", 0'; O2, OO 

-36 / Gt(0)Gt(0')Gt(0")-4i',2',3'^i.2.3 



ff >(0'i, 0'; 0, Oi)f r(02, 0"; 0', 02)f ^"'(Os, 0; 0", O 



if"=K;,+K-if3 

K'=K' +K' +K-K3-K2 



(A.4) 
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FIG. 2. Diagrammatic representation of the flow equation for the two-point vertex, see Eqs.(2.25) and (3.25). The left-hand 



side represents the derivative of the total irreducible vertex with respect to the flow parameter. The arrows represent the 
exact propagators, and an arrow with an extra slash represents the derivative of the cutoff function with respect to the flow 
parameter. 
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FIG. 3. Diagrammatic representation of the flow equation for the four-point vertex, see Eqs.(2.26) and ( i.28) 
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